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reach the current experimental upper bound, which can be tested by the experiments in the near 
future. 
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I. INTRODUCTION 



The Standard Model (SM) of particle physics, although greatly successful in phenomenol- 
ogy gives no explanations for parity (P) and CP violation. The observed neutrino oscilla- 
tions, the large baryonumber asymmetry and large energy density from non-baryonic dark 
matter (DM) in the universe are clear indications for newphysics beyond the SM. In the 
left-right (LR) symmetric models for weak interactions [l, 2|, [j, 0], the left- and right- 
handed fermions are treated equally, and P symmetry is restored prior to the spontaneous 
symmetry breaking (SSB). The LR models have other advantages such as the gauge group 
SU (2) L Cg> SU (2) R ® U (1)b-l can be elegantly embedded into grand unification theories, and 
the right-handed neutrinos are naturally required, etc.. 

The LR models may also contain DM candidates. In one of the minimal versions of the 
LR model which contains one Higgs bidoublet <fi and two Higgs triplets ^l,Ri a -^2 symmetry 
on the left triplet — > — A^ can be used to resolve the so-called VEV seesaw problem 
A direct consequence of this discrete symmetry is that the neutral component 5° L of Al 
can only annihilate or be produced by pairs, which makes 5° L a potential DM candidate. 
However, due to the fact that 5° L participates SU (2) L gauge interactions which is quite 
strong, the WIMP-nucleon elastic scattering experiments lead to a severe constraint on the 
dark matter relic density which is a few order of magnitudes below the observed value ij]. 
Thus it cannot be a main source of DM in the universe. 

A disadvantage of the minimal LR model is that although P can be broken spontaneously, 
the CP symmetry has to be broken explicitly, which looks quite unnatural. The reason is 
that without large fine-tuning on the Higgs potential the condition for spontaneous CP 
violation (SCPV) cannot be satisfied ilia 2|. Furthermore, in the minimal LR model 
with SCPV the predicted CP phase angle sin 2/3 ~ 0.1 in the unitarity triangle of the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix is far bellow the experimentally measured value 
of sin 2/3 = 0.671 ± 0.024 from the two B-factories jsj]. The minimal LR model also suffers 
from strong phenomenological constraints from low energy flavor-changing-neutral-current 
(FCNC) processes, especially the neutral kaon mixing which push the masses of the ri ght 



handed gauge bosons and some neutral Higgs bosons above the TeV scale 



Motivated by the requirement of both spontaneous P and CP violation, We have dis- 
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cussed an extension of the minimal LR model with two Higgs bidoublets [16j, [17|, [18|. In this 
two-Higgs-bidoublet LR model (2HBLR), the additional Higgs bidoublet may modify the 
Higgs potential such that the fine-tuning problem can be avoided. The extra Higgs bidoublet 
can also change the interferences among different contributions in the box-diagrams in the 
neutral meson mixin gs, and lower the bounds for right-handed gauge boson masses to be 



bellow the TeV scale 
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Note that the spontaneous P and CP violation in the LR models can also be useful for 
DM model-building. Before the SSB, the Lagrangian for the particle interactions prohibits 
the P-odd and CP-odd interactions, which may prevent the decays of the particles with odd 
CP parity. These particles can remain stable even after the SSB, provided that they do not 
develop VEVs and do not couple to the symmetry breaking sector. The simplest case would 
be that there is a gauge-singlet scalar field with odd CP parity, and has a vanishing VEV. 

In this work we discuss this possibility by considering an extension of the 2HBLR with a 
gauge-singlet complex field S which plays the role of DM candidate, and the stability of DM 



is purely protected by the discrete P and CP symmetries 



from the previous gauge-singlet models ( see, e.g. 19 



20 



21 



his model dis ting uishes itself 
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261 ] ) in that no 



ad hoc discrete symmetry of Z 2 type is introduced. This possibility has not been emphasized 
before in the literature, simply because most of the popular models such as SM and MSSM 
violate P and C maximally. This simple model shows that the DM may be connected to the 
fundamental symmetries of the quantum field theory. Recently, it is also noticed that the 
custodial symmetry of the gauge interaction can also be used to stabilize the DM candidate 



We calculate in this model the cross sections for the DM annihilation and the elastic 
scattering with the nucleus. The results show that for a large parameter space the DM relic 
density can be reproduced. The correlation between the DM annihilation and the elastic 
scattering off the nucleus depends on the Higgs and Yukawa sector of the model, and can be 
quite different from the ordinary gauge-singlet model which in some limit has a simple one 
to one correspondence. In general, the predictions for the direct detection experiments can 
be significantly larger and can even reach the current experimental upper bound for large 
Yukawa couplings. 

This paper is organized as follows: In Section. [Hj we outline the main feature of the 
model. In Section Hilt we discuss the parameter space, and give the formulas for main 
processes for the DM annihilation and the elastic scattering off the nucleus in a simplified 
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case where one Higgs bidoublet decouples from the theory. The case in which both Higgs 
bidoublet are active is discussed in Section [TV] We finally conclude in Section [V] 



II. THE LR SYMMETRIC MODEL WITH A GAUGE-SINGLET 

We begin with a LR model with a gauge group SU(2)l <S> SU(2)r <g> U{1)b-l which 
contains two Higgs bidoublets <fi (2,2*,0), x (2,2*, 0), a left (right) -handed Higgs triplet A L ( R ) 
(3(1), 1(3), 2), and a gauge-singlet 5(0,0,0) with the following flavor contents 




X 



Xi Xi 
X2 xl 



A 



L,R 



8+/V2 * 



L,R 



, S 



V2 



(So- + iS D ). 



(1) 



The introduction of Higgs bidoublet x is to overcome the problem of fine-tuning in generating 



the SCPV in 
constraints 



;he minimal LR model, and relax the severe low energy phenomenological 
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181 ] . Under the P and CP transformation, these fields transform as 





P 


CP 


X 


x' 


X* 


A L(R) 


A R(L) 


A* 


S 


s 


S* 



(2) 



We shall require P and CP invariance of the Lagrangian, which strongly restricts the struc- 
ture of the Higgs potential. For instance, for the terms involving the and A^^ fields the 



i 



most general potential is given by 

+A 1 [Tr(0V)]2 + A 2 {[Tr(( ^t 0)] 2 + [T r(^)} 2 } + A 3 [Tr(0V)Tr(00t)] 
+A 4 {Tr(0V)[Tr(0 f 0) + Tr(# f )]} 

+ Pl {[Tr(A L A[)] 2 + [Tr(A R A R )] 2 } + p 2 [Tr(A L A L )Tr(A[ A[ ) + Tr(A R A R )Tr(A R A R )] 
+p 3 [Tr(A L A L )Tr(A R A R )] + p 4 [Tr(A L A L )Tr(A R A R ) + Tr(A L A L )Tr(A R A R )] 
+ ai Tr(0t0)Tr(A L A[) + Tr(A R A R )] + a 2 Tr[(0V) + (#t)]Tr[(A L A L ) + (A R A R )] 
+a 3 [Tr(00 t A L A L ) + Tr(0VA R A R )] 

+/5![Tr(0A R 0tA[) + Tr(0tA L 0A R )] + /5 2 [Tr(0A R 0tA[) + Tr(0tA L 0A R )] 
+/3 3 [Tr(0A R 0tA[) + Tr(0tA L 0A R )], (3) 

where the coefficients p;, Aj, pi, on and $ in the potential are all real as all the terms are self- 
Hermitian. The Higgs potential V x a involving x field can be obtained by the replacement 
X «-> in Eq. ([3]). The mixing term V x ^a can be obtained by replacing one of by x i n ai l 
the possible ways in Eq. ([3]). 

In order to simplify the discussion, in this section we shall first consider a simple case in 
which the bidoublet x does not mix significantly with other fields. In this case the model 
is reduced to the minimal LR model plus a gauge- singlet, which already contains the main 
features of the complete model. We postpone the discussions on the x contributions into 
Section IIVI The most general Higgs potential involving the singlet field S is given by 

-V s = -j=p,&S + S*) - pl 2 s SS* - -^1{S + S*f + V2fx aS (S + S*)SS* 

+ -^-ft 3a (S + S*) 3 + \s(SS*) 2 - h*s(S + S*) 2 SS* - ±\ a (S + S*) 4 
6V2 4 lb 



-^=Ms + s*) + x i)S ss* - ±\ (s + s*) 2 



O l , (4) 



where 



1 = Tr(0t0), 2 = Tr(0t0 + 0t 0) and q 3 = T r(A[ A L + A^) . (5) 

Note that it only involves combinations of (S + S*) and SS*. The terms proportional to odd 
powers of (S — S*) are absent in the singlet self-interactions as they are P-even but C-odd. 



5 



Furthermore, they cannot mix with the Higgs multiplets in Oi because the three independent 
gauge-invariant combinations Oi(i = 1, ... ,3) in Eq. (jSJ) are both P- and C-even. Other 
possible Higgs multiplet combinations such as Tr(0^0 — 0^0) and Tr(A^A^ — A^A^j) are 
P-odd, thus cannot couple to S. The terms proportional to even powers of (S — S*) can be 
rewritten in terms of (S+S*) 2 and SS*. We have checked that the P and CP transformation 
rules for S defined Eq. ([2]) is acturally the only possible way for the implementation of the 
DM candidate. For future convenience we rewrite Vs in terms of the component field S a 
and So- 

— Vs = HoS a - -fJ^Sl - -^ 2 D S 2 D + -^oSl + HaDS a Sf) + -gX a S^ + -^\dS% + -\ aD S 2 a S 2 D 
+ j2(-^+^ S " + ^f S n)O i , (6) 
with the redefined coefficients 



1 r 



Hi,a — &i,<j, K,a — K,S ~~ \,D = K,S- (7) 

It follows from Eq. ([6]) that can only be produced by pairs, thus is a potential dark 
matter candidate. After the SSB, the Higgs multiplets obtain nonzero VEVs 

m = and <^> = ^ (8) 



where K\, k 2 , vl and vr are in general complex, and k = y | K il 2 + | ^2 1 2 ~ 246 GeV represents 
the electroweak symmetry breaking (EWSB) scale. The value of Vr sets the scale of LR 
symmetry breaking which is directly linked to the right-handed gauge boson masses. 

With the extra contributions from Vs to the whole Higgs potential V = V^a + Vs, one 
needs to redo the minimization with respect to <p and A^ ^. However, from Eq. (jHJ) and (J4j) 
it follows that the minimization conditions for and ^l(r) remain to have the same form 
as that in the minimal LR model. This is because the mixing introduced by the singlet S 
only changes the overall coefficients A*2 and ^ of the <p and A^ ; _r potential term in Eq. 
(J3J). Hence the mass matrix of the Higgs multiplet <fi and A^ r remains the same as that in 
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the minimal LR model in Refs. |28|, which also indicates that the additional potential 
term Vs in Eq. ([6]) does not help in resolving the fine-tuning problem. The fine-tuning can 
only be relaxed by introducing another bidoublet x- From the minimization condition for 
the singlets S a , one can eliminate one parameter \x a 

& = \y a - J2 (—m - XiAOi)) (9) 

where v a = (S a ) is the VEV of S a . In order to ensure the stability of the dark matter 
candidate Sd, we require that Sd does not obtain a nonzero VEV, (Sd) = 0, namely CP 
is not broken by the singlet fields. It follows that after the SSB, although P and CP are 
both broken, there is a Z2 type of discrete symmetry on Sd remaining in the gauge-singlet 
sector. The discrete symmetry is induced from the original CP symmetry. 

In the limit that ~ and k 2 <C K\ which comes from the phenomenology of neutrino 
masses and neutral meson mixings, the mass eigenstates for the Higgs bidoublet and triplets 
approximately coincide with the corresponding flavor eigenstates. The mass terms for the 
Higgs bosons and gauge bosons are listed in Table [H There is only one light SM-like Higgs 
h° from the real part of 0°, the mass of all the other scalars are set by vr which can be very 
heavy. 



Particles 


Mass 2 


Particles 


Mass 2 


h° = ffi 

Hi = ft 
A\ = 
Hi = 6% 
H$ = 5° L r 


m 2 = 2Xik 2 
m 2 H o = \cw? R 
m 2 A o = \a 3 v 2 R 
m 2 H o = 2 Pl v 2 R 

m H° = Up3~ 2 Pi) v r 


Hf = 0± 

= *r 
Ht = st 

H^ = 


m%± = ^a 3 v 2 R 

m] 1 ±± = 2p 2 v 2 R 

m 2 H ± = \(p 3 -2 Pl )v R 

m H±± = 5(03 -2pi)«| 
m 2 A o = \(p 3 -2 Pl )v R 


Zi 

z 2 


m li = m w 1 sec2 ®w 

2 _ g 2 v 2 R cos 2 8 w 
Z2 cos 29w 


wt = Wf 

wi = 


m% fri = g 2 K 2 /4 
m w 2 = 9 2 v 2 R /2 



TABLE I: The mass spectrum for the Higgs and gauge bosons in the LR symmetric model in the 
limit vl — and K2 <C k±. 4>® r and (ffi stand for real and imaginary component of 4>®. The gauge 
boson Z\(W\) corresponds to the Z(W) boson in the SM. 
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The mass terms for Sd and S a are given by 



M 2 D 




+ ^ + — + — <0*> , 


Ml 


= 2\ a vl + (jl 3a + 3jl a s) v a + ^-^ 


fa 




= «(-//l |(T + Xl^Va) , 




K? 


= 2k(-/2 2 , ct + A 2l0 -f CT ) , 






= v R (-p, S)(7 + Aa,^) , 





where M^o r , M^ r and M^ s0r denote the mixing between singlet and the other three 
neutral Higgs bosons. From the Lagrangian in Eq. (jSj) one can easily obtain the interaction 
terms among the scalars. Some of the relevant cubic and quartic scalar interactions are 
listed in Table [Til 



Interaction 


Coupling 


Interaction 


Coupling 


Interaction 


Coupling 


Interaction 


Coupling 


SdSdHH* 
S D S D h°H° 
S D S D AA* 


Al,D 

2A 2j d 

A3,D 


S D S D h° 
SdSdH® 

Q Q fjO 


Al,£)K 

+ Ai )CT k 
2A 2 ,dk 
2(-/i2, CT + A 2)CT k) 


SdSdSq- 

HH*S a 

hP 2*^ S a 
AA*S a 


2(m<tD + A ct du ct ) 

-/ii )(T + Ai^tv 

2(-^ 2 ,a + A 2i(7 U (J ) 

— M3,ct + A 3i(T t; ct 


SdSdH® 

c c uO 

J a ° 'a 12 2 

q Q Q 

h°h°H° 


h,DVR 

6X a v a 
aiv R 



TABLE II: The cubic and quartic scalar couplings between Higgs singlets and multiplets, where 
HH* stands for any states of {h°h°, H°H^ , A\A\ , H^H^ } and AA* stands for any states of 

{2*22*2, a" ao , Htnz,Hi+nz-,i$n°, h++h T }. 



III. DM IN THE LR SYMMETRIC MODEL 



There are a number of free parameters in the model such as the coefficients in the poten- 
tials and the VEV for S a . As shown in Eq. (flOi) . the mass of Sd is related to two energy 
scales v a and vr since (O3) ~ v 2 R . In the minimal LR model with the spontaneous CP vio- 



lation, the VEV vr of the right-handed 



from the K, B meson mixings 
interactions 30, 31 



9, 



11 



Hi ggs trip 



12 



13 



14 



et Ar is subjected to strong constraints 
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29j as well as low energy electroweak 



321 ] . The kaon mass difference and the indirect CP violation quantity 
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ex set a bound for vr around 10 TeV [30|, |33j, |34|. For a successful cold DM candidate Sd 



its mass should be roughly between 10 GeV and a few TeV for annihilation cross section of 
approximately weak strength. For simplicity here we consider a case in which v a is heavy 
~ Vr ~ 10 TeV, and before the SSB the Lagrangian has an approximate global U(l) 
symmetry on S, i.e. under S — > e lS S, which suppresses some of the parameters, namely 

Akb J^ui P-crSi /^3ct? fii,a ^ v ai 

Ao-S? A CT , Aj iCT <C 1, (11) 

which leads to a relatively light in comparison with vr, as it is the would-be Goldstone 
boson in the limit of exact global 17(1) symmetry. For a light Sd with Md < (9(lTeV), 
without significant fine-tuning, one needs £13^ /v a < 0.01 from Eq. (flQl) . It follows from Eq. 
( TlTT) that the coefficients for the quartic couplings and mixing terms in the potential V5 are 
roughly at the same order of magnitude 

A,j — — Kd, 

Ka ^ K,D. (12) 

One of the implications of this parameter region is that the mixing between S a and the 
SM-like Higgs h° will be small. This is because the mixing angle 9 is proportional to 

Thus the constraints from the precision electroweak data from LEP experiments becomes 
weak. There are of course other possible parameter regions. However, one will see in the next 
section that the parameter space corresponding to the approximate global U(l) symmetry 
leads to the correct magnitude of the relic dark matter density. 

A. Annihilation cross section 

The relic density of the gauge-singlet DM Sd can be calculated from the annihilation cross 
sections which depend largely on mass spectrum of the particles in the model, especially the 
mass of the DM candidate. Sd can be very lig ht. For 3 GeV < mo ^ 8 GeV which 
is consistent with the recent DAMA results 35[. Sd pairs can only annihilate to light 



fermion pairs through intermediate SM-like Higgs boson h°. In this case, there is an one to 




So 



Sn / 



h°,H° 




FIG. 1: Feynman diagrams for two DM candidate Sd annihilating into W\W\jZ\Z\ and tt final 
states. 



one correspondence between the DM relic density and the WIMP-nucleon elastic scattering 



cross section. The ratio between the two only depends on the mass of DM [22j, l36|. In order 
to satisfy both the DM relic density 0.105 < tt DM h 2 < 0.117 Q and the WIMP-nucleon 
elastic scattering cross section in the range 3 x 10 _41 cm 2 < o^ 1 <5x 10 _39 cm 2 reported by 
a large h Sp Sn coupling is inevitable, which may cause the invisible decay of 



DAMA 
h° produced at LHC 



36). 



Here we consider a different parameter rang in which Sd is heavier than the SM-like 
Higgs and in a mass range 200 < vtid < 500 GeV which can be covered by the CDMS and 
other experiments. Since we assume v a ~ vr ~ 10 TeV, most of the scalars are heavy except 
for the SM-like one. In this case, the possible annihilation products are h°h°, W\Wi/ 'Z\Z\ 
and fermion pairs qq, as shown in Fig. [1] and Fig. [2j The qq final states are dominated by 
the heavy t— quarks since the Yukawa coupling is the largest. For W\W\/ Z\Z\ finale states 
the only possible intermediate state is h°. For qq final states, the intermediate particles can 
be h° and H®. But Hf contribution is negligible as m H o ^> m^o. Since have nonzero 
B — L charge they can only couple to Majorana neutrinos. For a high Vr around a few TeV, 
the right-handed neutrinos are also heavy, which cannot appear in the final states. Thus 
in our model the dominant annihilation_processes in Fig. [1] are the same as in the minimal 



extension of SM with a gauge-singlet [22 
h°, H® 2 and 5*0- while the t— channel involves h° only. 

The relevant annihilation cross sections for Fig. [1] are given by 



26j | . For h°h° final states, the s— channel involves 



{4:EiE 2 (Jv) WlWl 
(AE 1 E 2 av) ZlZl 
{AE 1 E 2 av)a 



A? 



D 



m 



h° 



8tt 



m 



Wi 



+ 12- 



m 



Wi 



A? 



D 



1 



m 



h° \-2 



167T ' S 

3A?,p ml _ 
4tt s 



m 



12- 



m 



4m|. , i 



4m? 



|5fl 



m 



h° 



(14) 
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\ / 
\ / 

So W 



ft 



So \ / h° 

/ \ ft 



)•" — — -«( 



So/ 



h° 



5c I 

A D/ \ ft 

/ \ 



FIG. 2: Feynman diagrams for two DM candidate So annihilating into two SM-like Higgs bosons. 



and that for Fig. [2] is 

(4EiE 2 av) h o h o = 




+ 2A? D 



(m 2 D - m 2 h0 ) 2 

where s is the squared center-of-mass energy. E\ and E 2 are the energies of the in- 
cidental particles. The quantity Y is defined as Y(£) = arctanh(£)/£ with £ = 
y/(s — 4m^„)(s — 4m^)/(s — 2m^ — 2m 2 D ). For the cross section (4E 1 E 2 av) h o h o in Eq. 
f 1 1 5 p . the H^ffth scalar coupling a.\VR has been used 



a. 



B. Constraints from The DM Relic density 

The thermal-average of the annihilation cross section times the relative velocity (av) is a 
key quantity in the determination of the cosmic relic abundance of Sd- For nonrelativistic 
gases, (av) can be expanded in powers of relative velocity and x" 1 (x = mnjT). To the first 



order (av) ~ aox n , where n = 0(1) for s(p)-wave annihilation process 
formula for (av) is given by 39] 



3s 



(av) = a x~ 



m 



i) 



\2l)-u')x- 1 + 



The general 



(16) 



s/Am? D =l 



where u = E\E 2 av, and the prime denotes derivative with respect to s/Am 2 D . u and its 
derivative are all to be evaluated at s/Am 2 D = 1. The final DM density Q^Mh 2 is given by 



11 



38] 



Q DM h 2 = 1.07 x 10 9 f GeV- 1 (17) 

gj MpicT 

with 

x f = ln[0.038(n+l)(W^ 1/2 )^«^^o] 

-(n + 1/2) ln{ln[0.038(n + l)(g D / g l J 2 )M Pl m D a }} , (18) 

where Mpi = 1.22 x 10 19 GeV and gr> — 1 is the internal degree of freedom of Sp. is 
the total number of effectively relativistic degrees of freedom at the time of freeze-out. For 
particles playing the role of cold DM, the relevant freeze-out temperature is Xf = mp/Tf ~ 
25. Since we consider the range 200 GeV < m D < 500 GeV in our analysis, we obtain 
g* = 345/4. 




-0.20 -0.15 -0.10 -0.05 0.00 0.05 0.10 0.15 0.20 -0.20 -0.15 -0.10 -0.05 0.00 0.05 0.10 0.15 0.20 



FIG. 3: The allowed region of Xi^ and \r for different tud from DM relic density. The left panel 
corresponds to the annihilation involving only one Higss bidoublet. The right panel corresponds 
to the annihilation involving two Higgs bidoublet. See the text for detailed explanation. 

The total annihilation cross section uj is 

uj = (E^av)^ + (E 1 E 2 av) WlWl + (E l E 2 av) Zl z 1 + {EiE 2 av) tt ■ (19) 

From Eq. ffT5]) there are seven unknown parameters enter the expression of total annihilation 
cross section, namely, m h o, m^, \i,d, atiX^o, m^, m2 H o and vr. But u is highly insensitive 
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to m a and m 2 H o when s a and masses are around Vr. Furthermore the mass of H® can be 
related to Vr through m 2 H0 w 2piv R as it is shown in Table [H Thus only four parameters 

m h o, m D , A ljD and A/? = aiA 3)i j/(2pi) 

are relevant to our numerical analysis. In numerical calculations, we fix the mass of the 
SM-like Higgs to m^o = 120 GeV, and perform a numerical scan over the parameters Xr 
and A^d for the mass range 200 GeV < tud < 500 GeV. The allowed parameter space is 
shown in Fig. [3] (left panel), which gives an allowed range 

-0.18 < Ai D < 0.18 and - 0.48 < A^ < 0.47. 

The central region of this figure is excluded since these points cannot provide large enough 
annihilation cross section to give the desired relic abundance. For such a mass range of Sd, 
without significant fine-tuning, one needs A < 0.01 and jx/v a < 0.01 from Eq. (fTOj) . where 
A and fx denote the corresponding parameters in Eq. pip . Since the approximate global 
symmetry U{1) requires A <C Ai ; d, the region near Ai^ = in Fig. [3] is disfavored. 



C. Predictions for the DM direct detection experiments 



The current DM direct detection experiments, such as the CDMS 4Q] and XENON 4=11 ] . 
have imposed strong constraints on the WIMP-nuclen elastic scattering cross section for a 
wide range of DM mass. In our model, the DM candidate Sd interacts with nucleus M 
through Yukawa couplings interactions. For scalar interactions, the spin-independent elastic 
scattering cross section on a nucleus Af is given by 

4M 2 (AT) 



42 



43] 



7l 



(Zf p + (A-Z)f n 



(20) 



where M(N) = mpM^ j '{rriD + Mv) and Mj^ is the target nucleus mass. Z and A — Z are 
the numbers of protons and neutrons in the nucleus. f p ^ n is the coupling between WIMP 
and protons or neutrons, given by 

2 



fp,r 



f?q 

q=u,d,s 



(p.n) fTl 

a, 



p.n 



27 TG 



2^ <v 



q=c,b,t 



(21) 



where f$ = 0.020 ± 0.004, = 0.026 ± 0.005^/^ = 0.118 ± 0.062, /™ = 0.014 ± 0.003, 



?(p) 



(p) 



&] = 0.036 ± 0.008 and f^J = 0.118 ± 0.062 



(n) 



44j. The coupling between WIMP and 
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gluons from heavy quark loops is obtained from 

f<& = 1 - J2 ft ] ■ ( 22 ) 

q=u,d,s 

Traditionally, the results of WIMP-nucleus elastic scattering experiments are presented in 
the form of a normalized WIMP-nucleon scattering cross section a^ 1 in the spin-independent 
case, which is straightforward 

SI 1 M 2 (n) 

where Min) = moMnj (mn + M n ) is the reduced mass of the nucleon, and M n = m p , n 
denotes the nucleon mass. Because of f p ~ f n in our model 

AM 2 

7T 

The present bounds on the WIMP-nucleon elastic scattering cross section are a^ xp < 8 x 
1(T 44 cm 2 ~ 2 x 10~ 43 cm 2 from the CDMS [40J and a e n xp < 2 x 10~ 44 cm 2 ~ 4.3 x 10~ 43 cm 2 
for from the XENON [4l| for the DM mass range 200 ~ 500 GeV. 

The DM candidate Sjj interacts with nucleus M through their couplings with quarks 
by exchanging Higgs bosons h° and H®. Because H® is much heavier than h°, the main 
contribution comes from the /i°-exchange only. In this case, the WIMP-quark coupling a 



in Eq. (J2TJ) is given by 



a i = — 2 (25 ^ 

2v2m D m^o 

where (q = u, d, s, c, b, t) denotes the Yukawa coupling of the SM-like Higgs to the quarks 
with h^k/\^2 = m q . Using the allowed ranges for A^d we make predictions for the spin- 
independent WIMP-nucleon elastic scattering cross section cxf 7 . The numerical results are 
shown in Fig. 1 (left panel). One finds erf 1 < 7 x 10" 45 cm 2 for 200 GeV < m D < 500 GeV, 
which is far bellow the current experimental upper bounds. Nevertheless the future Su- 
perCDMS (Phase A) experiment 45J is able to cover part of the allowed parameter space, 
especially in the small mu region. 



IV. CONTRIBUTIONS FROM THE OTHER HIGGS-BIDOUBLET 

In this section, we generalize the previous discussions to the case in which the other 
bidoublet x mixes significantly with <p and A L>R . In this case the SCPV can be easily 

14 




FIG. 4: Predicted region of the spin-independent WIMP-nucleon elastic scattering cross section 
£T„ as functions of and Xi^d- The dashed line and solid line denote the present experimental 
upper bounds from the XENON and CDMS, respectively. The dotted line indicates the sensitivity 



of the future SuperCDMS [45|. The left panel corresponds to the annihilation involving only one 
Higss bidoublet. The right panel corresponds to the annihilation involving two Higgs bidoublet 
with the assumption of yq 1 ~ yq ~ . See the text for detailed explanation. 



realized 
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181 ] . Comparing with the previous case, the main differences are that 



there could be more scalar particles entering the DM annihilation and scattering processes. 
Furthermore, the new contributions from these particles may modify the correlation between 
the DM annihilation and WIMP-nucleus elastic scattering cross sections, which leads to 
significantly different predictions from the other gauge-singlet scalar DM models and the 
previous discussions. 

As shown in Eq. (1), the second Higgs bidoublet x contains two neutral Higgs particles 
Xi2- After the SSB, X12 m &y obtain VEVs iwi^/v^- The squared sum of all the VEVs 
including k 1)2 should still lead to v — (|^i| 2 + |^i| 2 + |u>i| 2 + IW2I 2 ) 1 / 2 ~ 246GeV. In the 
physical basis, some of the Higgs bosons from x could be light around electroweak scale. 
The number of the light Higgs depends on the Higgs potential. In most cases, there are two 
more light neutral Higgs rf( 2 and one pair of light charged Higgs rf 1 18| . This feature can 
be easily understand in the limit k-i ~ u>2 ~ 0. In this case, one can determine h°, rj® and 
7/2 from 0° an d Xi) an d V± from the mixing of (fif and xt- The number of kinematically 
allowed DM annihilation processes depends on the masses of the relevant particles. Here we 
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consider a case in which Sd is heavier than all the light scalars and the SM-like h° remains 
the lightest scalar, i.e. m^o < m v o ,m v + < mo- The quartic interaction and the s-channel 
annihilation in Fig. [2] now have seven possible final states. They are combinations of any 
two of the three neutral states 77^, 77°) and charged final states 7] + r]~. Note that each 
s-channel diagram in Fig. [2] may have h°, rj® and 77° as intermediate states besides S a and 
if 2- The t-channel diagram has six possible final states, due to the absence of the cubic 
scalar vertexes SdSdV^- 

The cubic coupling SpS^r}^ although can be different from that for S^Soh , may not 
modify the correlation between the DM annihilation and WIMP-nucleon elastic scattering 
cross sections in a significant way. As it is pointed out in Ref. [3^, the ratio R = (av)/a^ 1 
is highly insensitive to these couplings because they cancel out largely. R is only sensitive 
to the mass of DM candidate and the Yukawa couplings. In the minimal scalar DM model 
the Yukawa couplings are the same as that in the SM. it is shown that the value R scales as 
m 2 D [27]]. For small around a few GeV the value of R is in agreement with the DAM A 
results. A large mo around a few hundred GeV corresponds to a large R, which indicates 
that the WIMP-nucleon elastic scattering cross section may be far bellow the current direct 
detection bounds. 



In 2HBLR the Yukawa couplings can be significantly different from those in t 



re SM and 



the minimal LR model. Similar to the general two-Higgs-doublet model 



51 



52 



47 



49 



50, 



the Yukawa couplings are not simply determined by the quark masses. This is 
because with the introduction of the additional bidoublet, the fermion mass matrices and 
Yukawa matrices are not proportional to each other. In general the Yukawa couplings can 
be parameterized as y q = y/2£ q m q /v, the factor £ g depends on fermion flavor q and can 
be different from unity. For the DM annihilation processes, the heavy quark contribution 
dominates while for WIMP-nucleon scattering processes the light quarks are more important 
as the quark mass dependence are reduced in Eq. ff2T]) . For a large £ ? for light quark sector 
it is possible that the prediction for WIMP-nucleon elastic scattering cross section can be 
enhanced and the cross section for DM annihilation still coincides with the observed DM 
relic density. 

For a concrete numerical illustration, we choose all the masses m v o, m^o, m v ± at 180 GeV 
and keep m h o = 120 GeV. For cubic and quartic scalar couplings, we assume they are the 
same as that for the SM-like Higgs. Namely, the couplings of Sr>Sr>rji 2 an d SdSdVi 2V1 2 are 



1(3 



set equal to Xi t oV and A^d, respectively. Similarly, the cubic scalar couplings among the 
light Higgs particles h°, rfl 2 and V are se ^ equal to 3m^ /^. For a comparison we consider 
two sets of Yukawa couplings: 

i) All the couplings y^, yf, yf for h°qq, rj^qq, rj^qq respectively are nearly the same: 

yf ~ yf ~ yj. (26) 

with t/g = \/2m g /v. In this case annihilation cross section can be obtained simply by 
counting the number of new channels. In Fig. [3] (right panel), we give the constraints on 
for different m^. It is clear that there is a stronger constraint on the allowed 
parameter space, due to the increased number of intermediate and final states. For the 
WIMP-nucleon elastic scattering process, the WIMP-quark coupling a q in Eq. ff25l) is given 
by 

a Kdv ( yj | yf | y?\ 

q 2V2m D \m 2 h0 m 2 ^ m 2 o J 
Using the allowed Ai 5 £> and m,£> from Fig. [3l we calculate the spin-independent WIMP- 
nucleon elastic scattering cross section erf 7 . The numerical results are shown in Fig. |4] (right 
panel). We find that the erf 7 is enlarged roughly by a factor of three, which is however still 
below the current experiment upper bounds. 

ii) The Yukawa couplings for y^ 1 and y^ 2 are significantly larger in the light quark sector 
(for q = u,d, s) 

y m „ yf ~ 10yJ. (28) 

Since the annihilation process SdSd qq is dominated by heavy t-quarks, the enhanced 
Yukawa couplings yf and yf do not affect the total annihilation cross section. Thus the 
DM relic density remain unchanged. However, the predicted WIMP-nucleon scattering cross 
section erf 1 will be enhanced. The corresponding results have been shown in Fig. [5j We find 
that in this case erf 7 is enhanced by an order of magnitude compared with the one-Higgs 
bidoublet case. The future DM direct detection experiment SuperCDMS can cover most of 
the allowed parameter space. 

V. CONCLUSIONS 

In summary, we have discussed the possibility that the stability of DM can protected by 
the fundamental symmetries P and CP of quantum field theory. It can be realized in the 
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(27) 



200 250 300 350 400 450 500 

m D (GeV) 



FIG. 5: Predicted region of the spin-independent WIMP-nucleon elastic scattering cross section <t„ 
as functions of mo and \\,d in two Higgs bidoublet model with the assumption of yq 1 ~ yq 2 ~ lOy^ 
for three light quarks. 

framework of a generalized LR symmetric model which allows SCPV. The DM candidate in 
our model is a gauge-singlet which transforms under CP as an ordinary complex scalar. In 
this model no extra discrete symmetry is required. We have scanned the parameter space 
allowed by the relic DM density and made predictions for direct detection experiments. We 
found that the model was in agreement with the current measurement in a large parameter 
space. Based on the constrained parameter space, we have made predictions for the WIMP- 
nucleon spin-independent cross sections, and further studied the correlations with the DM 
annihilation. We have found that in this model the correlation could be significantly different 
from other gauge-singlet DM models. The DM-nucleon elastic scattering cross section could 
reach the current experimental upper bound for large Yukawa couplings for light quarks, 
which could be tested by the future experiments. 
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